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We put forward a dissipative preparation scheme for strongly correlated photon states. Our 
approach is based on a two-photon loss mechanism that is realised via a single four-level atom inside 
a bimodal optical cavity. Each elementary two-photon emission event removes one photon out of 
each of the two modes. The dark states of this loss mechanism are given by NOON states and 
arbitrary superpositions thereof. We find that the steady state of the two cavity modes exhibits 
entanglement and for certain parameters, a mixture of two coherent entangled states is produced. 
We discuss how the quantum correlations in the cavity modes and the output fields can be measured. 
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I. INTRODUCTION 

The performance of optical technologies such as 
metrology, communication and imaging can be improved 
beyond the limitations of classical physics if non-classical 
light sources are employed. The major challenge for the 
realisation of these quantum-enhanced schemes is the de- 
terministic generation of custom-tailored photon states 
for specific applications. For example, quantum infor- 
mation schemes based on continuous variable entangle- 
ment [l[ have the advantage that entangled light fields 
can be generated unconditionally, but high-quality re- 
sources with a large degree of entanglement arc diffi- 
cult to produce. Another example is given by optical 
intcrferomctry, a technique that is employed in applica- 
tions like gravitational wave detectors, laser gyroscopes 
or optical imaging. It usually aims at the precise esti- 
mation of a relative phase acquired by the light on its 
way through the interferometer. The achievable preci- 
sion of the phase estimation with classical light sources is 
bound by the standard quantum limit (SQL) and scales 
as 1/y/W, where N is the (mean) number of photons 
in the input field @]. On the contrary, a better preci- 
sion with the same amount of resources can be obtained 
with entangled light fields. A Prominent example for 
photon states that allow one to beat the SQL is given 
by so-called NOON states d i], which also give rise to 
phase super-resolution Realistic scenarios that in- 
clude photon losses of the interferometer require states 
with a more complicated structure than NOON states 
for breaking the SQL 0, • Another example is given 
by coherent entangled states @, (CES) that allow 
one to beat the SQL in lossless interferometers [TTI. |l2T[. 
These states yield a better phase estimation than NOON 
states in lossy interferometers fl3| if the same mean num- 
ber of input photons are taken into account. However, 
it remains challenging to produce NOON [TJ, [l5| and 
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states with a high success rate and with 
a large number of photons. 

One of the most successful techniques for the prepa- 
ration of quantum mechanical systems in a desired state 
is dissipation. Prominent examples are given by laser- 



and evaporative cooling that allow one to realize a Bosc- 
Einstein condensate. Recently, the concept of dissipa- 
tive quantum state preparation |18l - [20| was transferred to 
the many-body domain where dissipation alone prepares 
strongly correlated states. The challenge in dissipative 
quantum state preparation is to design a suitable dissipa- 
tive process £r such that the desired state \4>) is station- 
ary with respect to £r, he., Cr(\'4')('4'\) = 0- For exam- 
ple, a dissipative contact interaction was investigated in 
one-dimensional molecular [2ll - |24| and polariton [25l - [27j 
systems, where dissipation effectively results in a repul- 
sion between particles. The entanglement of two distant 
atomic ensembles via spontaneous emission was investi- 
gated in [2^, [2{|, and the simulation of open quantum 
systems with ion systems was considered in [30l. l3l|. 

Here we present a dissipative preparation scheme for 
strongly correlated photon states inside an optical cav- 
ity with two modes a and b. We engineer a two-photon 
loss term via a single, laser-driven four-level atom that 
couples to the cavity modes, see Fig. [1] Each elementary 
emission event induced by this two-photon loss term re- 
moves one photon out of mode a and one photon out of 
mode b. The dark states of the engineered dissipator £r 
are given by all strongly entangled NOON states and su- 
perpositions thereof, and the stationary state inside the 
cavity can be well approximated by a mixture of two CES 
states for specific parameters. We show that the steady 
state of the cavity modes alone is entangled, and point 
out how the entanglement of the cavity modes and of the 
output field can be measured. 

This paper is organised as follows. We give a detailed 
description of our system (see Fig. [T]) in Sec. |TTJ Here 
we also derive an effective master equation for the cavity 
modes alone, which we obtain by an adiabatic elimina- 
tion of the atomic degrees of freedom. In Sec. lIIII we anal- 
yse the steady state of this effective master equation and 
identify the dark states of the engineered two-photon loss 
mechanism. The entanglement of the two cavity modes 
is discussed in Sec. |IVI We employ the Negativity [32[ 
and an inequality [33| based on Einstein-Podolsky-Rosen 
type observables as sufficient entanglement criteria. The 
predictions of both measures are compared, and Sec . IVl 
indicates how the criterion based on the inequality [33| 
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FIG. 1. (Color online) A single four-level atom interacts with 
two cavity modes and a classical laser field. A- ln and Bi n 
are coherent input fields that drive the corresponding cavity 
resonantly, and A out , B out are the output fields. The inset 
shows the atomic level scheme. The laser field with frequency 
ul and Rabi frequency £Il couples to the 2) |3) transition, 
and the cavity mode with frequency ui a interacts with the 
|3) o |1) transition. The second cavity mode with frequency 
tjf, interacts with the |4) o |2) transition. g a and gt are the 
single-photon Rabi frequencies corresponding to mode a and 
b, respectively. The parameters jij are the decay rates of the 
various transitions, S and A label the detuning of the cavity 
fields with transition |3) o |1) and |4) -B- |2), respectively, 
and e is the two-photon detuning. 

could be measured experimentally. The latter Sec. |V] 
is mostly concerned with the entanglement of the out- 
put field. We identify two suitable modes of the output 
field whose entanglement can be inferred from two-mode 
squeezing spectra. Finally, we address the experimental 
realisation of our scheme in Sec. IVII and conclude with a 
summary and outlook of our results in Sec. IVII1 



II. REDUCED MASTER EQUATION FOR THE 
CAVITY MODES 

We consider a single four-level atom that interacts with 
two cavity modes and a classical laser field, see Fig. U It 
was shown theoretically [34j and experimentally [35j that 
the considered level scheme can give rise to a strongly 
enhanced Kerr nonlinearity. In addition, this level con- 
figuration allows one to engineer a two-photon absorp- 



tion process (2(| H(| . Note that we chose the two-cavity 
setup in Fig. [T] because it allows us to present our model 
in a clear and unambiguous way. However, our setup 
could very well be realised with a single cavity where the 
two modes can be either two different polarisation or fre- 
quency modes. In particular, the two modes could have 
the same frequency and orthogonal polarisations. 

The aim of this section is to derive an equation of mo- 
tion for the reduced density operator g-p of the two cav- 
ity modes. We begin with a detailed description of the 
system shown in Fig.[TJ The input field A- m (-Bin) is a co- 
herent field that resonantly drives the cavity mode with 
frequency cu a (u>b). The corresponding Hamiltonian is 

H m = m* a e iu ^a + hQ,* b e lU)ht b + H.c., (1) 

where H.c. stands for the Hermitian conjugate, and a 
(b) is the annihilation operator of the cavity mode with 
frequency u a (ujb). The Rabi frequencies fl a and fib are 
determined by the input power of the fields A m and B m , 
respectively. The cavity mode with frequency u> a couples 
to the atomic transition |3) <-> |1), and the mode with 
frequency ujb interacts with the atom on the |4) <-> |2) 
transition. In rotating- wave approximation (RWA), the 
interaction of the atom with the cavity modes is described 
by the Hamiltonian 

H c = -% a a|3)(l| - % 6 6|4}(2| +H.c, (2) 

where g a (gb) is the single-photon Rabi frequency on the 
1 3} <-> |1) (|3) O |2)) transition. The detuning of the first 
cavity mode with the |3) <-> |1) transition is denoted by 
5, and A is the detuning of the second mode with the 
1 4) <-> 1 2) transition, 

5 = uj a - oj 31 , A = uj b - uj 42 . (3) 

The resonance frequencies on the |3) o |1) and |4) O |2) 
transitions have been labelled by W31 and W42, respec- 
tively. In addition, the atom interacts with a classical 
laser field with frequency wl and Rabi frequency JIl, and 
this field couples to the |3) 4-> |2) transition. In rotating- 
wave approximation, the atom-laser interaction reads 

H L = -M7 L |3)(2|e-' iWLt +H.c. (4) 

The free time evolution of the cavity modes and of the 
atomic degrees of freedom is given by Hp and Ha, re- 
spectively, 

+ huj b b%, (5) 
H A = R[wa|2>(2| +w 3 |3>{3| +W4|4>(4|], (6) 

and we set uj\ — in Eq. With these definitions, we 
arrive at the master equation for the combined system of 
the atomic degrees of freedom and the two cavity modes, 

q= - l -\H F + H in + Hp L + H^ + H c , e \+C 1 Q + C K Q. (7) 
n 
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The term C^g in Eq. ([7]| describes spontaneous emission 
of the atom and is given by 
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(8) 



where jij is the full decay rate on the transition \i) <-> \ j) 
(see Fig. [T]). The atomic transition operators are defined 



S+ = |3><1|, 5+ = |3)(2|, 5 3 + = |4)(2|, (9) 

and S~ = {S^~y . The last term C K g in Eq. (J7]) accounts 
for photon losses at the cavity mirrors and reads 



C K g = -^-(a) ag + go) a — 2aga)) 

(tfbg+ gtfb-2bgtf), 



Kb 



(10) 



where n a (/sj) is the damping rate of mode a (b). 

In appendix|21 we derive from Eq. (UJ) the master equa- 
tion for the reduced density operator g F of the cavity 
modes alone, 



Of = Tr A g = Qn + g22 + g33 + 044 , 



(11) 



and g vv denotes If the two-photon detuning e = 

uj a — ojl — ^2 vanishes, the master equation for the density 
operator gp of the cavity modes in an interaction picture 
with respect to Hp in Eq. (0) is given by 



£?f = Ci a QF + C K g F + C v g F + £rg~F, 



where 



(12) 



(13) 



accounts for the coherent driving of the cavity modes 
and the cavity decay term C K gF is defined in Eq. (|10p . 
The two terms CuSf and CrgF in Eq. (|12j) represent the 
influence of the single atom on the evolution of the two 
cavity modes. More specifically, 



£uQf = —i U[{ba)^ba, 



Qf\ 



(14) 



describes a coherent two-particle interaction between 
photons in modes a and b. The strength of this inter- 
action is determined via the parameter 



U 



\9a\ 2 \g b 
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.-A. 



(15) 



Note that the sign of U depends on the sign of the detun- 
ing A defined in Eq. It follows that the interaction 
described by Cjj can be repulsive or attractive. The last 
term in Eq. (|12[) is given by 

C-f&f = —77 \(ba)%agY + QF(ba)'ba — 2bag-p(ba) 1< ~\ (16) 



and represents a two-photon loss term with decay rate 



r = 



\9aW 



(A2+ 7 f 2 /4)|fi L p 



742- 



(17) 



The dissipator Cr gives rise to the emission of correlated 
photon pairs. In each elementary emission process, Cr 
removes one photon out of mode a and one photon out of 
mode b. In the following, we focus on the quantum corre- 
lations between the cavity modes that are induced by the 
two-photon loss term Ct- For the rest of this paper, we 
will thus consider the purely dissipative scenario where 
the conservative photon-photon interaction vanishes, i.e. 
U = 0. This situation arises when mode b is resonant 
with the 1 4) o |2) transition and thus A = 0. 

Next we summarise the conditions for the validity of 
the reduced master equation (|12j) . The adiabatic elimi- 
nation of the atomic degrees of freedom requires that the 
atomic decay rates 7^ are large as compared to the other 
system parameters, 



742,731,732 > K a ,K b , y/N a N b T, o , Q b , \6\, (18) 

where N a (N b ) is the largest relevant photon number in 
mode a (b). The master equation (|12[) for the cavity 
modes holds under the assumption that the two-photon 
detuning e vanishes. Since e = can be adjusted only 
within a certain accuracy, we establish conditions that 
ensure the validity of Eq. (|12p for e 7^ in Appendix [B] 
Finally, the restriction to terms up to third order in the 
expansion of the density operator (|A7|) requires 



N a \g a \ 



< 1 



N b \g b 



<1, 



(19) 



where N a (N b ) is the maximal photon number in mode 
a (b). 



III. STEADY STATE ANALYSIS 

In this Section we characterise the steady state of 
the master equation (|12l) with a purely dissipative two- 
photon interaction (A = U = 0). Throughout this work 
we assume that the cavity decay rates n a and n b are dif- 
ferent from zero. In this case and for time-independent 
input fields £l a and tt b , the master equation (|12[) exhibits 
a unique steady state. If the two-photon decay term T 
vanishes, this unique steady state is the pure state [l8| 
g F = IV'XV'I, where \ip) = \a a ,a b ) and \a a ) (\a b )) de- 
notes a coherent state in mode a (b) with amplitude 
a a = -2iQ a /n a (a b = -2iQ. b /K b ). 

Next we turn to the general situation with r / 0. 
In this case, the analytical steady state of the master 
equation (fT2"j) is difficult to obtain. A numerical study 
shows that the steady state of Eq. (fT2"|) is in general a 
mixed state. In addition, an intuitive understanding of 
the structure of this state can be gained via the dark 
states [lH \D) of the two-photon loss term in Eq. ([To]) 
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FIG. 2. (Color online) Population of the cavity Fock states \n a ,nt) in steady state for (a) o = f2j, = 1.32 x 10 2 742, 
n a = n b = 10- 2 742, (b) n a = ft b = 6.26 x 1CT 3 742, K a = K b = 0.5 x 10- 2 742, and (c) Q a = fi 6 = 1.16 x lO"^, 
k„ = Kt = 10 _3 742. The common parameters in (a)-(c) are V = 742/IOO, A = U = 0, 5 = and e = 0. The mean photon 
number for all states in (a)-(c) is (N) = 5. 



that obey Ct{\D){D\) = 0. Since Cr removes one photon 
out of each cavity mode in every elementary emission 
event, all pure dark states \D) must be of the form 



\D) 



Cl|0a,O) +C2|0,^6) 



(20) 



where \<p a ) and \ipb) are arbitrary states (including the 
vacuum) of mode a and b, respectively. Note that Cr 
supports an infinite number of dark states, and the most 
general dark state is given by a mixture of different pure 
dark states that are of the form of \D) in Eq. (|2T)|) . In 
the regime where the two-photon decay rate T dominates 
over the coherent drive terms fl a , fib and the cavity de- 
cay rates Ki, K2, one can expect that the stationary state 
of Eq. (fT2"f is approximately a dark state of Cr alone. 
This is confirmed by Fig. [5] that shows the population of 
the cavity Fock states \n a , n&) in steady state for various 
parameters. These results were obtained via a numeri- 
cal solution of Eq. (fT2")l . The relative importance of the 
two-photon loss term (fl~6| is the smallest in Fig. [Ha) and 
the largest in Fig. HJc). In the latter case it is apparent 
that the steady state is comprised of the dark states in 
Eq. (|2U)) . Only Fock states |0, rib) and \n a ,0) are signif- 
icantly populated, while the population of other states 
\n a ,rib) with n a ^ and rib 7^ is strongly suppressed. 
This result can be understood as follows. All Fock states 
\n a ,rib) with n a ^ and rib 7^ experience not only 
the cavity loss term C K , but the additional (strong) two- 
particle losses Cr- This suppresses the population of 
these states via the cavity input fields. 

The cavity pump fields fl a , Of, induce transitions be- 
tween dark states |n a ,0), |0, rib) and neighbouring Fock 
states \n a , 1), |l,nt} (n a ^ and rib 0)- The latter 
states are not dark states of Cr and thus decay rapidly. 
However, we point out that this mechanism does not in- 
duce an indirect decay of the dark states |0, rib), \ n a, 0) if 
r becomes much larger than all other system parameters. 
The reason is that the transitions |0, rib) ~ * \^-i n b) — > 
\Q,n b - 1) (|n ? ,0) -> \n a ,l) -> \n a - 1,0)) occur at an 
effective rate fl 2 a /T < 1 (fig/r < 1) for fl a < T 
(fib <C T) and are therefore negligible. This result can 
be regarded as a manifestation of the quantum Zeno ef- 



fect [22| and explains the sharp population contrast be- 
tween dark states \n a ,0), |0,n&) and neighbouring Fock 
states \n a ,l), \l,n b ) in Fig.^c). 

The most general pure dark state in Eq. (f2"0)) is a coher- 
ent superposition of the states |</> a ,0) and |0, tfb)- This 
is a remarkable feature since it implies that all strongly 
entangled NOON states 



1 



|NOON) = ^=(\N,0) 



|0,iV» 



(21) 



and coherent superpositions thereof are dark states of the 
dissipator Cr- Next we investigate whether the steady 
state of Eq. (fT2"]) contains entangled dark states or just a 
mixture of the states |</> a ,0) and \0,ipb)- Note that this 
information is not contained in Fig. [2] because it does not 
contain any information about the off-diagonal matrix 
elements of the density operator. We find that the steady 
state exhibits non- vanishing coherences in the Fock basis 
if the modes a and b appear in a completely symmetric 
fashion in the master equation (|12[) . For example, for 
the parameters in Fig. [2][c), the diagonalisation of the 
numerically computed density operator shows that the 
steady state is very well approximated by a mixture of 
two coherent entangled states (CES), 

gr w Pl |CES + (ai))(CES + (ai)| 

+ Pa|CES_(a 2 )><CES_(a 2 )|, (22) 

where px w 0.499, p 2 « 0.459, a x -2.23i, a 2 ~ -2.29i 
and the CES states are defined as [1, [l3| 



|CES±(a)) 



v/2(l ± e~l«l") 



(|a,0)±|0,a)). (23) 



The overlap between the approximate state in Eq. 
and the full numerical solution in the trace norm is 
95.44%. 

Recently, the performance of CES in interferometric 
precision measurements was discussed in fl3| . Since it 
was found that these states perform better than NOON 
states in a lossy interferometer, we discuss the suitabil- 
ity of the state in Eq. (f2"2")l for interferometric precision 
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measurements. We find that the mixture in Eq. (f2"2")l 
performs approximately as good as classical light with a 
well-defined phase. Note, however, that the performance 
of a pure CES and the mixture in Eq. (f2"2"j) differ only if 
the losses of the interferometer are smaller than approxi- 
mately 25%. In this regime of a near-perfect interferom- 
eter, pure CES entangled states allow one to access the 
quantum regime that is currently inaccessible with our 
preparation method. 

The previous discussion indicates that the steady state 
of the two cavity modes can be prepared in an entangled 
state. The entanglement properties of the cavity modes 
and of the output field are discussed in the following Sec- 
tions |IV] and El respectively. 

IV. ENTANGLEMENT OF THE CAVITY FIELD 

The two modes of the cavity form a bipartite quan- 
tum system. By definition, the quantum state gp of the 
cavity field is said to be entangled if and only if it is 
non-separable, and gp is separable if and only if it can 
be written as 

&F = S ^PiQa(j)®Qb{3)- (24) 

3 

Here Q a (j) and Qb{j) are normalised states of the modes 
a and 6, respectively, and the parameters pj > are 
constrained by JN Pj = 1 ■ 

We employ two criteria that are both sufficient for the 
entanglement of the cavity modes. We begin with the 
Negativity [32| that is defined as 

M(QF) = \(l-\\e^\\i), (25) 

where g^ a denotes the partial transpose of gp with re- 
spect to the subsystem a. The trace norm ||A||i = 
Tr(vMtA) of an operator A is equal to the sum of its 
singular values. Note that the Negativity is an entan- 
glement monotone, and hence j\f(gp) > is a sufficient 
criterion for the entanglement of the two field modes. 

The second criterion is derived in [33| and states that 
the system is in an entangled quantum state if the to- 
tal variance of two Einstcin-Podolsky-Roscn (EPR) type 
operators u and v of the two modes satisfy the inequality 

f + (^cav) 2 ) < 2 , (26) 

where 

Ucav = X a + X b , Vcav = p a - p b . (27) 

For any operator X, we define SX = X — (X). and Xk 
and pk are local operators which correspond to mode k € 
{a, 6} with frequency ujk- The only restriction imposed 
on the operators Xk and pk is that they must obey the 
commutation relation 

[xk,Pi] = iSki , (28) 
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FIG. 3. (Color online) Numerical evaluation of (a) the Neg- 
ativity and (b) the inequality in Eq. ()33|) for the two cavity 
modes in steady state. The results are shown as a function 
of the two-particle loss rate V and for different intensities of 
the coherent driving fields. The parameters are <f> = 7r/2, 
K a = K b = 1(T 2 742, A = U = 0, S = and e = 0. The red 
dashed line corresponds to 51 a = fit, = 0.2 x 10~ 2 742, the blue 
dotted line stands for il a = SI5 = 0.4 x 10 _2 742, the green 
dashed-dotted line corresponds to Q, a = Q b — 1 x 10 _2 742 
and the black solid line stand for Q, a = Q b — 1.6 x 10 _2 742- 

but are otherwise arbitrary. Here we choose for ik and 
pi the quadrature operators of the cavity fields 

s„ = 4=(5e~^ + 5 t e < *) > (29) 
v2 

p a =^Xae- i ^-~a\e i % (30) 

5 6 = 4=(6e~ < * + 6 t e^), (31) 
v 2 

p b = JL(&e-^-&+e^), (32) 

where a = e IWat a and b = e ZUJbt b are slowly varying in 
time. Since u cav and u cav can be identified with the op- 
erators corresponding to the centre-of-mass motion and 
the relative momentum of two quantum mechanical oscil- 
lators, respectively, they are called EPR type operators. 
The phase <j> is arbitrary and will be optimised such that 
a maximal violation of the inequality (|26[) is achieved. 
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The inequality in Eq. (f2"6]l is equivalent to 

(: [(^cav) 2 + (^>cav) 2 ] :> < , 



(33) 



where : : denotes normal ordering (all creation operators 
to the left) with respect to the operators a, b and their 
adjoints. For reasons that will become clear later (see 
Sec. EJ, we will employ Eq. ([3"31 instead of Eq. 
With the help of Eqs. $Q and P9" ]) -([g2" l> . we find 



[(a cav ) 2 + (<5f) cav ) 2 ] :) =2[{5tf5a) 



(34) 



We numerically solve for the steady state of the den- 
sity operator of the two cavity modes via the master 
equation (|12p for various parameters T and different in- 
tensities of the coherent input fields. The result for the 
entanglement criteria in Eqs. (|25|) and ([34]) is shown in 
Figs.^a) and (b), respectively. It follows that the steady 
state of the system exhibits entanglement. The Negativ- 
ity and the normally ordered variance of the EPR-typc 
operators show a qualitatively similar (but mirrored) be- 
haviour in the weak driving regime. A different situation 
arises if the Rabi frequencies f2 a and f2& become com- 
parable to the cavity decay rates. Although the cavity 
modes are entangled for all parameters in Fig. [3ta) (the 
Negativity is larger than zero), the sufficient entangle- 
ment criterion in Eq. (|33[) is not fulfilled for larger values 
of r/742. We show in Sec. IVl and Appendix [Cl that the 
quantity in Eq. ((34)) and hence the entanglement crite- 
rion in Eq. (|33p can be measured experimentally. The 
discussion above shows that this approach allows one to 
capture the entanglement for a large range of parameters, 
but it fails to detect the entanglement in some cases. 



V. ENTANGLEMENT OF THE OUTPUT FIELD 



The output field of the cavity is a multi-mode field. 
The question whether the output field is entangled thus 
requires to specify the corresponding modes. In standard 
input-output theory [H, |39[ , the output fields A out and 
Bout [see Fig. Q] are defined as 



^out(t) = 
B ou t(t) = 



1 



2 71 

1 



/2tt 



dcj e -M*-ti),4(a;,ti), (35) 



Here A(uj,ti) and B(u,ti) are the Heisenberg operators 
of the continuous output modes taken at time t\ — > 00. 
The latter operators obey the equal-time commutation 
relations 

[A(uj,t),A^uj',t)) = [B{oj,t),B\u\t)\ = S(u-<J). 

(37) 

Here we focus on the entanglement between two modes 
corresponding to the central frequencies uj a and LOb of 



A out and B out , respectively. Our aim is to follow a simi- 
lar approach as in Eq. (|3"3"|) where we employ the variance 
of EPR-type operators as a sufficient criterion for entan- 
glement [33J . Therefore, we have to define position- and 
momentum-like operators that obey the canonical com- 
mutation relation (|28[), In order to achieve this, we have 
to construct a discrete mode of A out from the continuous 
mode operator A(uj,ti). This transition can be achieved 
if we average A{u, t\) over a small frequency interval Alo 
centered at u> ai 



Aui/2 



An = 




V2?rAw 



dwe'^+^MK+w.ti), (38) 



dtu / d^I out (t)e i 



-Alj/2 -00 



where the output field A out (t) = e , " a *^4 out (t) is slowly 
varying in time. Similarly, we define a discrete mode of 
the output field i? Q ut centered at 



Bo 



1 



Aw/2 



dcj / dtB out (ty 



(39) 



-Acu/2 



and B ont {t) = e %u)bt B out (t). The modes An and £> rep- 
resent modes of the output fields A out and i? ut that are 
experimentally accessible via spectral filtering. Further- 
more, we note that (AlAa) can be interpreted [itj as 
the mean number of photons emitted into the mode An 
within the time 1/Aw (the same statement holds for the 
mode Bo). 

Since the two discrete modes An and Bq obey the com- 
mutation relations 

[A ,Al] = [B ,Bl] = l, [A o ,Bl] = [Ao,B o } = 0, (40) 

we can define position- and momentum-like operators for 
each mode that obey the canonical commutation rela- 
tion 



1 
1 



x A = ^={A Q e-^+J^e^), 



PA = ^.{A Q e-^-Ale^) 
x B = ±J,B Q e- i * + Ble i <t>), 



(41) 
(42) 
(43) 
(44) 



The phase in Eqs. (|4"Tj) and (j4"4"|) is arbitrary and later 
on chosen such that the necessary separability condition 
is maximally violated. We can now proceed as in Sec. IIVI 
and define EPR-typc operators 



Uo = x A + X B , 



Vn=PA-PB- 



(45) 
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-4-2 2 4 

lj/k 



FIG. 4. (Color online) Squeezing spectra. The blue dashed 
line corresponds to S u (w), the red dotted line to S v (uj) and 
the black solid line shows S u (w) + S v (uj). The parameters are 

(f> = 7T/2, K = K a = Kb = 1CT 2 742, T = 742/IOO, O a = fib = 

0.7 x 1CT 2 742, A = (7 = 0, <5 = and e = 0. 



With the help of Eqs. (gSJ) and (|H ]) -(|3g ]) , we find that 
the normally ordered total variance of the operators uq 
and vq is given by 



Auj/2 




du[S u (uj) + S v (uj)], 



-Aw/2 



(46) 

where S u (u>) and S v (uj) are related to the two-mode 
squeezing spectra of the output fields. A definition of 
these functions and how they can be calculated numer- 
ically and measured experimentally is provided in Ap - 
pendix [Cl According to the entanglement criterion in [33( , 
the modes Ao and Bq are entangled if the left-hand side 
of Eq. (|4"rJ)) is smaller than zero. If Aw is sufficiently 
small, this means that modes Ao and Bq are entan- 
gled if S u (0) + S v (0) < 0. The numerical evaluation 
of the squeezing spectra S u and S v in Fig. [4] demon- 
strates clearly the entanglement of the output modes. 
For the chosen parameters, we find a two-mode squeez- 
ing of — 1.3dB in S u and — 6.1dB in S v at the central 
frequency ui = 0. 

So far we considered the entanglement between the 
two central frequency components at u> a and u>b- Note 
that our approach can be generalised to modes centered 
around u> a — u and uj^ + lj if the output spectra of modes 
a and b are completely symmetric and if the two modes 
are interchangeable in the master equation (|12j) . In pre- 
vious works P, lU |42[ , other EPR-type operators than 
those in Eqs. (g5]) with (gTJ-flUJl were defined, and their 
variances can also be related to squeezing spectra. How- 
ever, in the latter approach it is not obvious to identify 
the modes of the output field that are entangled. 

Finally, we address the experimental verification of the 



entanglement inside the cavity. To this end, we note that 

oo 

- / du,[S«(u) + S v (lj)} = (: [(5u cav ) 2 + (5v cliv ) 2 } :), 



and hence it follows that the field inside the cavity is en- 
tangled if the sum of S u (oj) + S v (uj), integrated over all 
frequencies, is negative [see Eq. ([55)) ]. Note that the vari- 
ance of the cavity fields can be measured directly without 
recording the squeezing spectra [43| . 



VI. EXPERIMENTAL REALIZATION 

Next we discuss the experimental implementation of 
our scheme in Fig. [TJ A key requirement of our scheme is 
that the single atom should have a very well defined po- 
sition that changes very little over the range of an optical 
wavelength. Ideal candidates are therefore single trapped 
neutral atoms [HI, |4f| or ions j46H50l | inside an optical 
cavity. Recently, several experiments with 40 Ca + ions in- 
side a high-finesse optical cavity have been reported (46l — 
l49j . A suitable candidate would also be given by a 138 Ba + 
ion where the level scheme in Fig. Q] could be realised 
between two J = 1/2 Zecman manifolds via polarisation- 
and frequency selection. Alternatively, the level scheme 
in Fig. Q] can be realised with artificial atoms [f3l[ coupled 
to microwave fields. 

In our approach we obtain an effective, engineered mas- 
ter equation for the cavity modes alone via an adiabatic 
elimination of the atomic degrees of freedom. In partic- 
ular, this requires that that the cavity decay rates are 
much smaller than the atomic spontaneous emission de- 
cay rates, see Eq. ([T5)l . In a recent experiment (4(| H3| 
with a 40 Ca + ion in a high-finesse cavity, the ratio be- 
tween the atomic decay rate and the corresponding cav- 
ity decay rate is j/k m 16. At the same time, the cou- 
pling constant g was larger than the atomic decay rate, 
g/j ~ 2 [4rjj . However, we show in appendix [D1 that 
it is advantageous to realise our scheme with g/j < 1. 
This opens up the possibility to increase the length L of 
the cavity which will reduce the magnitudes of the cou- 
pling constant \g\ ~ 1/y/L and of the cavity decay rate 
k ~ 1/L. From the above example, we conclude that 
values of j/k ps 100 should be achievable with current 
technology. Valid choices for the remaining system pa- 
rameters that comply with the conditions (|18[) and (fTI)]) 
of our model are discussed in appendix [Dj In summary, 
an increase in the number of photons in the cavity re- 
quires a reduction in the two-photon decay rate T. Note 
that this does not necessarily result in a smaller entangle- 
ment between the cavity modes. On the contrary, Fig. [3] 
indicates that the maximal entanglement is attained for 
smaller two-photon decay rates T if the number of pho- 
tons increases (see Sec. II VI) . 
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VII. SUMMARY AND OUTLOOK 

This paper investigates a scheme for dissipative quan- 
tum state preparation for two optical cavity modes a and 
b. We engineer an effective reservoir via a single, laser- 
driven four level atom that interacts with both cavity 
modes on separate transitions. The adiabatic elimina- 
tion of the atomic degrees of freedom gives rise to a mas- 
ter equation for the cavity modes alone and contains a 
two-photon loss term. Each elementary emission event 
induced by this two-photon loss term removes one pho- 
ton out of mode a and one photon out of mode b. We 
find that the dark states of this loss term are given by 
all entangled NOON states and superpositions thereof. 
If the two cavity modes are interchangeable in the corre- 
sponding master equation, the steady state of the cavity 
field exhibits entanglement. We employ the Negativity 
as well as an inequality [HI based on Einstein-Podolsky- 
Rosen type observables as sufficient entanglement crite- 
ria. While the former is an entanglement monotone, the 
latter can be measured via two balanced homodyne de- 
tection setups of the output fields. Furthermore, we de- 
fine suitable modes of the output fields and show that 
they can be entangled as well. The entanglement of the 
output modes can be verified experimentally via the mea- 
surement of two-mode squeezing spectra. 

The stationary state inside the cavity can be well ap- 
proximated by a mixture of two CES states for specific 
parameters. While a pure CES state is able to access 
the quantum regime in intcrferomctric precision measure- 
ments with small photon losses inside the interferometer, 
the mixture prepared by our dissipative scheme does not 
perform better than classical light with a well defined 
phase. However, modifications of our scheme may pave 
the way towards the dissipative preparation of NOON 
states and CES that are relevant for quantum-enhanced 
technologies. In the ideal case, it may help to prepare 
photon states that allow one to achieve the ultimate 
quantum limit Q in interferometric precision measure- 
ments with realistic photon losses. 
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Appendix A: Derivation of the reduced master 
equation 

Here we outline the derivation of the effective master 
equation for the cavity modes in Eq. (|12p . The starting 
point for our derivation is the full master equation in 
Eq. (0. Next we apply a unitary transformation W = 
W F ® Wa to Eq. ©, where W F = cxp[iH F t/h] acts only 



on the cavity modes, and 

W A =exp[i(H A A + e|2)(2| +5|3)(3| + (A + e)|4)(4|)t] 

(Al) 

acts only on the atomic degrees of freedom. The density 
operator in the new frame is denoted by g = WgW^ and 
obeys the equation of motion 

Q= - l -[H + H c ,~o\ + C-,~q + Cio~q, (A2) 

where 

H = -%|2>(2|+<y|3)(3| + (A + e)|4)(4| 

+ n L |3)(2|+n£|2)(3|], (A3) 

and e = Lu a — cjl — ^2 is the two-photon detuning. The 
super-operator C\o = An + £ K hi Eq- (|A2[) accounts for 
the external driving via the input fields and the damping 
of the cavity modes, where £; n and C K are defined in 
Eqs. (TIB")) and (fTU|) . respectively The master equation for 
the transformed density operator g-p of the cavity modes 
is obtained if we trace over the atomic degrees of freedom 
in Eq. (lA"2ll . 

Qf = (ig* a [a\ hi] + i9b[&, Q42} + H.c.) + & Qf- (A4) 

In order to eliminate the coherences g^i an d f?42 from 
Eq. 1|A4[) . we solve Eq. (|A2|) pcrturbativcly in the Hamil- 
tonian He describing the interaction between the atom 
and the cavity modes. In order to obtain the desired 
expansion of the full density operator g in the coupling 
constants g a and gi,, we re- write Eq. (|A2|) in a form where 
the atom-cavity interaction is separated from the other 
terms, 

1 

q = Cqq- -[H c ,q], (A5) 

and the super-operator Co is defined by 

1 

Cog = --[H ,g} + C 1 g + dog. (A6) 
Expansion of the density operator in Eq. (|A5|) as 

e = f> (fe) , (A7) 

fc=0 

where g^ denotes the contribution to g in kth order 
in He, leads to the following set of coupled differential 
equations 

^ (0) = Cog^ , (A8) 
^ =C a ~g^ - l -[H c ,~g (k - 1] ], k>0. (A9) 

Equation (|A8|) describes the interaction of the atom with 
the classical laser fields to all orders and in the ab- 
sence of the cavity fields. Higher-order contributions to 
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g can be obtained if Eq. (|A9|) is solved iterativcly. Equa- 
tions (|A8j) and (|A9|) must be solved under the constraints 
Tr(£(°>) = 1 and Tr(f>W) = (k > 0). We employ 
a Markov-type approximation and assume that the atom 
reaches its steady state on a timescale that is fast as com- 
pared to the typical evolution time of gp induced by the 
atom-cavity coupling. In addition, we suppose that the 
cavity decay rates K a , Kb and the Rabi frequencies fi a , fib 
are small as compared to the atomic decay rates 7y and 
hence neglect the contribution of the super-operator Cio 
to Co- Under these conditions, the set of equations (|A8j) 
and (|A9|) can be solved in a straightforward manner if 
Co is represented by a matrix. However, the procedure 
is tedious since He contains the operators a and b, and 
hence one has to keep track of the operator ordering. 
Up to third order, g w + fjW + g^ + g( 3 \ and for 
vanishing two-photon detuning e = we find 



031 



Ag a \gb\ 2 b^bag F , 



042 = Ag b \g a \ 2 bag F a\ 



where 



A 



(A10) 
(All) 



(A12) 



(A + ' J742 /2)|tt L | 2 - 

Finally, we substitute Eqs. (|A10|) and (jAllI) in Eq. jSH) 
and obtain the master equation (|12p. 



Appendix B: Finite two-photon detuning £ / 

Here we establish conditions that ensure the validity of 
Eq. (|12[) for non-zero two photon detuning. For simplicity 
we limit the discussion to the case of a purely dissipative 
photon-photon interaction (U = A = 0). In the ideal 
case e = 0, only the third-order term contributes to the 
coherences in Eqs. (|A10|) and (|A11[) . The first-order term 
for e ^ gives rise to an additional term C\g F on the 
right-hand side of Eq. flT2l . Up to second order in e, we 
find 

£i£»f = - r[Hi, Qf] 
n 



(a'agp + gpa'a — 2agpa i ), (Bl) 



where 



#1 = -%o| ^TT^m a T « 



ri=73 pr7 £ ' 



(B2) 
(B3) 



and 73 = 731+732 is the full decay rate of state |3). Here 
Hi is just a frequency shift of mode a, and the term pro- 
portional to Ti describes an additional decay channel for 
photons in mode a. This term can be neglected provided 
that 



Fl=7 W <<KQ - 



(B4) 



The second-order contribution g^ is equal to zero for 
all parameters. On the other hand, other third-order 
terms occur for £ ^ that give rise to other two-particle 
processes in Eq. (fT2"|) . The magnitude of these terms 
relative to the terms proportional to T in Eq. (fT2"j) is 
negligible provided that |e|, \5\ <§C 731,732,742 and 



|e|(73 +742) 



< 1, 



|fe| 



< 1. 



(B5) 



Appendix C: Definition and measurement of 
squeezing spectra 

The two-mode squeezing spectra are defined as 

00 

S u (u>)=2 / dr cos ujt(T : 5u out (t + T)Su out (t) :), (CI) 



S v (oj)=2 j drcoswr(r : Sv ont (t + T)Sv out (t) :), (C2) 


where 

"out = x A + x Bl v out =pa~Pb, (C3) 
and the quadratures of the output field are given by 



PA = -^r(i O ute- i0 - Al ut e 



(C4) 
(C5) 
(C6) 



p B = ^{B out e-^ - Bl^). (C7) 

The operator T in Eq. (|C2|) orders products of annihi- 
lation operators such that their time arguments increase 
from right to left, and products of creation operators are 
ordered such that time arguments increase from left to 
right. Since the definition of S u and S v comprises only 
expectation values of normally and time-ordered prod- 
ucts of the output fields, we can easily represent S u and 
S v in terms of the cavity fields [38j ] . 



S u (uj) = 2k j drcoso;r(T : 5u ca , v (t + r)(5u cav (i) :), 





(C8) 



S v (u>) — 2k J dTCOswr(T : Sv cav (t + r)(5w cav (t) :), 

(C9) 

where u cav and v cav are defined in Eq. (j27]) . The general 
setup for the measurement of the two-mode squeezing 
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FIG. 5. (Color online) Measurement of the two-mode squeez- 
ing spectrum in Eq. ()C8[) . The output field A ovA (B out ) 
is superimposed with a strong coherent field at a 50 : 50 
beamsplitter, and the photocurrents of each detector are sub- 
tracted. The photocurrents of the two balanced homodyne 
detections are added and fed into a spectrum analyser. The 
resulting spectrum is given by Eq. (|C10[) and directly propor- 
tional to S u (u)- 



spectra [4l], |42j is shown in Fig. [5j The output field A out 
(-Bout) is superimposed with a strong coherent field at 
a 50 : 50 beamsplitter, and the photocurrents of each 
detector are subtracted. A similar measurement is per- 
formed for the output field -E> ou t- The local oscillator 



fields are given by A\ oc — \a\e and B 



where 



|a| = |/3 1 and the phase 9 — 4> — tt/2 determines the phase 
4> occurring in the definition (|C5j) - (|C7|) of the quadrature 
operators. The photocurrents of the two balanced homo- 
dyne detections are added and the resulting current i is 
fed into a spectrum analyser. The recorded spectrum is 
then given by [52| 



P u (oj) =- I dr coswt lim [{i(t)i(t + r)) - {i{t)f 

7T J t— >O0 




= -P s hot[l + r}S u (ijj)], 



(C10) 



where P s hot is the shot-noise limit and 77 is the detection 
efficiency. The spectrum S v (w) can be measured with 
the same setup if the phases of the local oscillators are 
changed according to 8 ^ 9 + n/2, and the photocurrents 
of the two homodyne detections have to be subtracted 



rather than added. 



Appendix D: Parameter choices 

The validity of our approach requires that the condi- 
tions in Eqs. ([Tg]) and (flU]) are fulfilled. Here we dis- 
cuss the possible choices for the system parameters that 
comply with these conditions. For simplicity, we as- 
sume that the coupling constants, cavity decay rates, 
Rabi frequencies and the Fock state cutoffs are iden- 
tical for both modes. We thus set \g a \ = \gb\ = \g\, 
K\ = «2 = K, O a = Oj, = O and N a = Nf, = N . Further- 
more, all atomic decay rates are assumed to be the same, 
7 = 731 = 732 = 742, and we introduce dimensionless 
parameters p = p/7 that are denoted by a tilde. Since 
we are only interested in the purely dissipative situation, 
we limit the discussion to the case A = such that the 
conservative photon-photon interaction U vanishes. We 
begin with the cutoff N in the Fock state basis that is 
determined by the Rabi frequency £1 driving the cavity 
modes, the cavity decay rate re and the two-photon decay 
rate T. Note that Eq. (|18[) requires these parameters to 
be much smaller than unity, and thus N cannot be ar- 
bitrarily large for realistic values of k. Next we discuss 
the condition for the two-photon decay rate V in Eq. (fT5|) 
and the inequality in Eq. (|T9l) . For a given cutoff N of 
the Fock state basis, these two conditions read 



AT< 1, 



< 1. 



(Dl) 



The two small parameters x and NT determine the abso- 
lute values of the coupling constant g and of the Rabi fre- 
quency corresponding to the external laser that drives 
the ion. With the explicit expression for T in Eq. (|17l) . 
we find 




I fir 



ivvr 

2x 



(D2) 



Since x and NT are both small parameters, one can 
achieve \g\ < 1 and therefore the strong coupling regime 
of cavity QED is not required. This is advantageous be- 
cause Eq. (|D2|) implies that |Ol/<7| = y N/x, and hence 
|Ol| can be much larger than \g\. It follows that small 
values of \g\ ensure that the scaled Rabi frequency I^lI 
remains reasonably small. In addition, a smaller coupling 
constant facilitates the realisation of cavity decay rates 
that are much smaller than the atomic decay rates (see 
Sec. EH- 
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